Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120-km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.     

Eocene tectono-thermal rejuvenation of an upper Paleozoic-lower Mesozoic terrane in the Cordillera de Carabaya, Puno, southeastern Peru, revealed by K---Ar and Ar/Ar dating

K---Ar dates for muscovites and biotites in granitoid rocks and hydrothermal ore deposits of the northeastern parts of the plutons making up the Triassic Carabaya batholith, underlying the axial Cordillera Oriental of northern Puno Department, southeastern Peru, are markedly variable and mutually discordant. Steep transverse gradients are defined in the apparent ages of both micas, which decrease systematically from SW to NE, delimiting a ca. 25-km-wide, longitudinal zone of anomalously young Mesozoic to Paleocene dates. Age minima of 37±1 Ma are attained in three of the four studied transects. ⁴⁰Ar/³⁹Ar step-heating analyses of selected micas confirm the occurrence of a thermal disturbance, and modeling of the spectra suggests that argon loss in muscovites attains at least ca. 75% in the northeastern part of the zone. A single K-feldspar spectrum yielded a minimum at 31 Ma, and apatite fission-track age cluster at ca. 31 and 18.5 Ma. The affected granitoid rocks generally display little megascopic evidence of tectonism, but microscopic deformational fabrics increase in intensity with apparent decreasing K---Ar age, paralleling a marked increase in alkali feldspar ordering. Secondary fluid inclusions trapped within the microfabrics reveal that the plutonic rocks were penetrated by a homogeneous H₂O---CO₂---CH₄---NaCl fluid at ca. 300–400°C and 0.7–2 kbar. This fluid is implicated in the degassing of the rocks. These diverse data are interpreted as evidence for a major, but moderate-temperature (400°C) and brief, tectono-thermal event at ca. 37±1 Ma (biotite closure temperature)—i.e., at the Eocene-Oligocene boundary. The K-feldspar ⁴⁰Ar/³⁹Ar data and the Oligocene fission-track dates may record the later stages in the event, whereas the Miocene fission-track dates are tentatively ascribed to a distinct Neogene episode. Essentially identical geochronological and petrological relationships have been documented in the Cordillera Real of northwestern Bolivia by McBride et al. (1987), permitting the delimination of a disturbed belt paralleling the South American plate boundary and more than 450 km long. The tectono-thermal domain, which we term the Zongo-San Gabán Zone, constituted the foreland boundary of the Andean orogen in the vicinity of the Arica Deflection during the late Eocene Incaic orogeny. This regional thermal event, which involved the basement, appears to have resulted from compressional or, in some segments, transpressional tectonics.     

Time scale of hydrothermal water-rock reactions in Yellowstone National Park based on radium isotopes and radon

We have measured ²²⁴Ra (3.4 d), ²²⁸Ra (5.7 yr), and ²²⁶Ra (1620 yr) and chloride in hot spring waters from the Norris-Mammoth Corridor, Yellowstone National Park. Two characteristic cold-water components mix with the primary hydrothermal water: one for the travertine-depositing waters related to the Mammoth Hot Springs and the other for the sinter-depositing Norris Geyser Basin springs. The Mammoth Hot Springs water is a mixture of the primary hydrothermal fluid with meteoric waters flowing through the Madison Limestone, as shown by the systematic decrease of the (²²⁸Ra/²²⁶Ra) activity ratio proceeding northward. The Norris Geyser Basin springs are mixtures of primary hydrothermal water with different amounts of cold meteoric water with no modification of the primary hydrothermal (²²⁸Ra/²²⁶Ra) activity ratio. Using a solution and recoil model for radium isotope supply to the primary hydrothermal water, a mean water-rock reaction time prior to expansion at 350°C and supply to the surface is 540 years assuming that 250 g of water are involved in the release of the radium from one gram of rock. The maximum reaction time allowed by our model is 1150 years.     

Andesite and dacite genesis via contrasting processes: the geology and geochemistry of El Valle Volcano,Panama

The easternmost stratovolcano along the Central American arc is El Valle volcano, Panama. Several andesitic and dacitic lava flows, which range in age 5–10 Ma, are termed the old group. After a long period of quiescence (approximately 3.4 Ma), volcanic activity resumed approximately 1.55 Ma with the emplacement of dacitic domes and the deposition of dacitic pyroclastic flows 0.9–0.2 Ma. These are referred to as the young group. All of the samples analyzed are calc-alkaline andesites and dacites. The mineralogy of the two groups is distinct; two pyroxenes occur in the old-group rocks but are commonly absent in the young group. In contrast, amphibole has been found only in the young-group samples. Several disequilibrium features have been observed in the minerals (e.g., oscillatory zoning within clinopyroxenes). These disequilibrium textures appear to be more prevalent among the old- as compared with the young-group samples and are most likely the result of magma-mixing, assimilation, and/or polybaric crystallization. Mass-balance fractionation models for major and trace elements were successful in relating samples from the old group but failed to show a relationship among the young-group rocks or between the old- and young-group volcanics. We believe that the old-group volcanics were derived through differentiation processes from basaltic magmas generated within the mantlewedge. The young group, however, does not appear to be related to more primitive magmas by differentiation. The young-group samples cannot be related by fractionation including realistic amounts of amphibole. Distinctive geochemical features of the young group, including La/Yb ratios〉15, Yb〈1, Sr/Y〉150, and Y〈6, suggest that these rocks were derived from the partial melting of the subducted lithosphere. These characteristics can be explained by the partial melting of a source with residual garnet and amphibole. Dacitic material with the geochemical characteristics of subducted-lithosphere melting is generated apparently only where relatively hot crust is subducted, based on recent work. The young dacite-genesis at El Valle volcano is related to the subduction of relatively hot lithosphere.     

Lead isotopic evidence for evolutionary changes in magma-crust interaction,Central Andes,southern Peru

Lead isotopic measurements were made on Andean igneous rocks of Jurassic to Recent age in Moquegua and Tacna Departments, southernmost Peru, to clarify the petrogenesis of the rocks and, in particular, to investigate the effect of crustal thickness on rock composition. This location in the Cordillera Occidental is ideal for such a study because the ca. 2 Ga Precambrian basement rocks (Arequipa massif) have a distinct Pb isotopic signature which is an excellent tracer of crustal interaction, and because geomorphological research has shown that the continental crust was here thickened drastically in the later Tertiary.Seven samples of quartz diorites and granodiorites from the Ilo and Toquepala intrusive complexes, and seven samples of Toquepala Group subaerial volcanics were analyzed for Pb isotopic compositions. The plutonic rocks range in age from Jurassic to Eocene; the volcanic rocks are all Late Cretaceous to Eocene. With one exception, the Pb isotopic ratios are in the ranges ²⁰⁶Pb/²⁰⁴Pb= 18.52–18.75, ²⁰⁷Pb/²⁰⁴Pb= 15.58–15.65, and ²⁰⁸Pb/²⁰⁴Pb= 38.53–38.74. The data reflect very little or no interaction with old continental material of the Arequipa massif type.Lead from four Miocene Huaylillas Formation ash-flow tuffs, two Pliocene Capillune Formation andesites and five Quaternary Barroso Group andesites has lower ²⁰⁶Pb/²⁰⁴Pb than that in the pre-Miocene rocks, but relatively high ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb (²⁰⁶Pb/²⁰⁴Pb= 18.16–18.30, ²⁰⁷Pb/²⁰⁴Pb= 15.55–15.63, ²⁰⁸Pb/²⁰⁴Pb= 38.45–38.90). Tilton and Barreiro [9] have shown that contamination by Arequipa massif granulites can explain the isotopic composition of the Barosso Group lavas, and the new data demonstrate that this effect is evident, to varying degrees, in all the analysed Neogene volcanic rocks. The initial incorporation of such basement material into the magma coincided with the Early Miocene uplift of this segment of the Cordillera Occidental [32], and thus with the creation of a thick crustal root. The data strongly imply a relationship between crustal thickness and degree of crustal contamination of magmas in this area, but a rigorous relationship is not yet established.     

The latitudinal distribution of carbon monoxide across the pacific from California to Antarctica

The results of a research study of the carbon monoxide concentration from California to 90° S, Antarctica are presented. The data both extend and support other research studies of the latitudinal distribution of carbon monoxide in that higher concentrations are evident over the Northern Hemisphere than over the Southern Hemisphere. Carbon monoxide concentrations range between 50 to 60 ppb with a few peaks into the 60s in the latitudinal area south of the ITCZ and values of 80 ppb or higher at latitudes north of Hawaii. A comparison is also made of carbon monoxide and ozone concentrations along the flight tract between California and Antarctica, over the Ellsworth Mountains of Antarctica, and between 78° S and the South Pole. These ozone-carbon monoxide data show statistically significant negative correlations in the upper troposphere and lower stratosphere over Antarctica. It is believed that this is a good indication of mixing across the tropopause.